Approximate quantum error correction can lead to better codes 
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We present relaxed criteria for quantum error correction 
wliicii are useful wlien tiie specific dominant quantum noise 
process is known. These criteria iiave no classical analogue. 
As an example, we provide a four-bit code which corrects for a 
single amplitude damping error. This code violates the usual 
Hamming bound calculated for a Pauli description of the error 
process, and does not fit into the GF(4) classification. 



I. INTRODUCTION 

Quantum error correction is the reversal of part of the 
adverse changes due to an undesired and unavoidable 
quantum process. Code criteria for perfect error correc- 
tion have been developed [0-|4j, and all presently known 
codes satisfy these criteria exactly. Furthermore, most 
of these codes belong to a general GF(4) classification 
1^. In these schemes, quantum bit (qubit) errors are de- 
scribed using a Pauli (/, X, Y, Z) error basis, and coding 
is performed to allow correction of arbitrary unknown er- 
rors. 

However, in the usual case in the laboratory, one works 
with a specific apparatus with a particular dominant 
quantum noise process. The GF(4) codes do not take 
advantage of this specific knowledge, and may thus be 
sub-optimal, in terms of transmission rate and code com- 
plexity. Unfortunately, there is no general method to 
construct quantum codes in a non-Pauli basis, and few 
such codes are known |^-[ic|] . The code criteria are gener- 
ally very difficult to satisfy, and without the Pauli basis, 
no way is known to apply classical coding techniques for 
quantum error correction. 

In this paper, we develop a novel approach to quan- 
tum error correction based on approximate satisfaction 
of the existing quantum error correction criteria. These 
approximate criteria are simpler and less restrictive, and 
in certain cases, such as for amplitude damping, codes 
can be found relatively more easily. The approximate 



criteria also admit more codes, therefore, codes requiring 
shorter block lengths may also be possible. 

For example, using this approach we have discovered 
a four bit code which corrects for single qubit amplitude 
damping ]lO| , pd| errors. Such a short non-degenerate code 
is impossible using the Pauli basis. The reason is that the 
effects 13| to be corrected in a non-degenerate code 
have to map the codeword space to orthogonal spaces if 
the syndrome is to be detected unambiguously. Hence, 
the minimum allowable size for the encoding space is the 
product of the dimension of the codeword space and the 
number of effects to be corrected. The single qubit am- 
plitude damping effect operators expressed in the Pauli 
basis are: 



(I + v/1-7)/ + (1 - v/l-7)f^a 



V7 



(1) 
(2) 



These describe the changes due to the loss of zero or 
one excitation to the environment. To first order in the 
scattering probability 7, n-\-l possible effects may happen 
to an n-qubit code using the Aq, Ai error basis, so it 
follows that n > 3 is required. In contrast, in the Pauli 
basis, any ax or ay error must be corrected by the code, so 
that 2n-\-l possible effects must be dealt with. It follows 
that at least n > 5 is required for a non-degenerate Pauli 
basis code, in contrast to the four bit code which we 
demonstrate in this paper. 

The lessons are that (1) better codes may be found 
for specific error processes, and (2) approximate error 
correction simplifies code construction and admits more 
codes. Approximate error correction is a property with 
no classical analogue, because it makes use of slight non- 
orthogonalities possible only between quantum states. 
We describe our approach to this problem by first ex- 
hibiting our four-bit example code in detail. We then 
generalize our results to provide new, relaxed error cor- 
rection criteria and specific procedures for decoding and 
recovery. We conclude by discussing possible extensions 
to our work. 
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II. FOUR BIT AMPLITUDE DAMPING CODE 

Consider the single qubit quantum noise process de- 
fined by 



E{p) = ^^P^ 
1 







V7 




(3) 
(4) 



known as amplitude damping. The probabihty of losing 
a photon, 7, is assumed to be small. To correct errors 
induced by this process, we encode one qubit using four, 
with the logical states 



1 



|0^) = _ 10000) + 11111) 



1 
71 



10011) + |1100) 



(5) 
(6) 



A circuit for encoding the logical state is shown in Figj^. 
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FIG. 1. Circuit for encoding a qubit. The third mode con- 
tains the input qubit. The rotation gate in the first qubit 
performs exp(— i7r(Tj;/4). 

The possible outcomes after amplitude damping may 
be written as 



(7) 



where k are strings of O's and I's serving as indices for 
each error, and we use the notation Aoio -- = AqA\Aq ■ ■ ■. 



[■) is convenient shorthand |11| for a mixed state (tensor 
sum of un- normalized pure states). In the following, the 
squares of the norm of |-) states will give their proba- 
bilities for occurring in a mixture. For the input qubit 
state 



IV'...) =a|0L)+6|li 



(8) 



all possible final states occurring with probabilities 0(7) 
or above are 



\(f>oooo) = a 
+ b 

I'/'iooo) = 



10000) + (1-7)2|1111) 



V2 

(1-7)[|0011) 



(9) 



11100)] 



7(1 - 7) 
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a(l-7)|0111) + 6|0100) 



|0oioo) 
|0ooio) 
l^oooi) 



/7(1- 


7) 
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/7(1- 


7) 
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'7(1- 


7) 



a(l-7)|1011) + 6|1000) 
a(l-7)|1101) + 6|0001) 
a(l-7)|1110)-H6|0010) 



The usual criteria used to study quantum codes (re- 
viewed in section HI) require that (Oi|^|A^|OL) = 

(iLlAUfelU), but 



(Ol|AJooo^oooo|Ol) = 1 - 27 -I- 37' -f 0(7') 
(1l|4ooo^oooo|1l) =1-27 + 72. 



(10) 
(11) 



So the code we have constructed does not satisfy the 
usual criteria. We will demonstrate that the code sat- 
isfies new approximate error correction conditions later 
on, and revisit the recovery procedure afterwards. First, 
we exhibit how it works. 



A. Decoding and Recovery Circuit 

Let us denote each of the four qubits by \ni),. . .,|n4). 
Error correction is performed by distinguishing the five 
possible outcomes of Eq.(^), and then applying the ap- 
propriate correction procedure. The first step is syn- 
drome calculation, which may be done using the circuit 
shown in Fig. ^A. There are three possible measurement 
results from the two meters: (M2, M4) = (0,0), (1,0) and 
(0,1). Conditioned on (M2, M4), recovery processes Wk 
implemented by the other three circuits of Fig. || can be 
applied to the output \n1n3) from the syndrome calcula- 
tion circuit. 
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FIG. 2. (A) Circuit for error syndrome detection. The 
measurement result is used to select out of three actions. 
If the result (M2, M4) is 00, 10, or 01, circuits (B), (C), or (D) 
are applied, respectively, to recover the state. The angles 9, 9' 
are given by tan 9 = {1 — 7)^ and cos = 1 — 7. The rotation 
gate and controUed-rotation gate specified by the angle 9 per- 
form the functions exp{i9ay) and Ai{exp{i9ay)) respectively 
in the notation of |l4|. 
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If (M2, Mi) = (0,0), then jning) is: 
|00> + (1-7)^111)1 , J(l-7)(|01) + |10)) 



72 



V2 



(12) 



To regenerate the original qubit, the circuit of Fig. gB is 
used: a controUed-not is apphed using I713) as the control, 
giving 



a|0) 



+ (1-7)^11) 
V2 



6|1) 



(1-7)(|1) + |0)) 



V2 



(13) 



|ni) can now be used as a control to rotate jna) to be 
parallel to |0). We obtain as the final output: 



(1 j^i 



|0) (14) 



(l-7)(a|0)+6|l)) + 0(7^)10) |0> , (15) 



with the corrected and decoded qubit left in |7ii) as de- 
sired. 

If (Af2, Mi) = (1,0), the inferred state before syndrome 
measurement is (/>iooo or (/)oioo- In either case, \nin^) is 
in a product state and the third qubit has the distorted 
state: 



1^3) 



(1 _ ^)^ 



a(l-7)|l) + fo|0) 



(16) 



To undo the distortion, we apply the non-unitary trans- 
formation in Fig. The combined operation on jria) [n-^] 
due to the NOT gate, the controUed-rotation gate and the 
measurement of the ancilla bit can be expressed in the 
operator sum representation: M{p) = NopN^ + NipNl, 
where 



iVo = |0)(1| + (1 - 7)|1>(0|, iVi = V7(2-7)|l)(0| ■ 

(17) 

The Nq and iVi operators correspond to measuring the 
ancilla to be in the |0) and |1) states respectively. If the 
ancilla state is |0), we obtain the state: 



(1 _ ^)3^ 



a|0) +6|1) 



(18) 



in the third mode because Nq preferentially damps out 
the &|0) component in Eq.(p^). We get an error message 
if the ancilla is in the |1) state. Finally, if {M2, M^) — 
(0,1) the same procedure can be applied as in the (M2, 
M4) — (1,0) case, with ni and swapped. 

The fidelity, defined as the worst (over all input states) 
possible overlap between the original qubit and the recov- 
ered qubit is 



1-57^+0(7'), 



(1-7)37' 



(19) 
(20) 



Note that the final state Eq.(15) is slightly distorted. 
This occurs because the recovery operation is not ex- 
act, due to the failure to satisfy the code criteria exactly. 
Furthermore, the circuits in Fig. and ^ have a finite 
probability for failure. However, these are second order 
problems, and do not detract from the desired fidelity 
order. 



III. APPROXIMATE SUFFICIENT CONDITIONS 

We now explain why our code works despite its vio- 
lation of the usual error correction criteria. The reason 
is simple: small deviations from the criteria are allowed 
as long as they do not detract from the desired fidelity 
order. This section presents a simple generalization of 
the usual error correction criteria which makes this idea 
mathematically concrete. These approximate error cor- 
rection criteria are sufficient to do approximate error cor- 
rection. We expect that a more complete theory giving 
necessary and sufficient conditions to do approximate er- 
ror correction is possible, but have not obtained such a 
theory. Nevertheless we hope that the simple sufficient 
conditions presented here will inspire other researchers to 
develop a general theory of approximate error correction. 

Quantum error correction is performed by encoding 
logical basis states in a subspace C of the total Hilbert 
space H so that some effects of £ can be reversed on 
C. Let the noise process be described in some operator 
sum representation £{p) = J^neK ^npAl^, where /C is the 
index set of A, the set of all effects An appearing in the 
sum. We denote by Are C A the reversible subset on C, 
and let /C^e = {^^| A^ € Are} be the index set of ^re- 
in other words, the process £'{p) = X^neK ^nP-^n is 
reversible on C. Are includes all the effects satisfying the 
condition (l]-||] 



(21) 



where Pc is the projector onto C, and gran are entries of 
a positive matrix. When Eq.(|2|) is true for some subset 
of effects An which form an operator-sum representation 
for the operation £, there exists some operator-sum rep- 
resentation of £ such that £' {p) = X^neyc- ^npA^^ and: 



PcAlAnPc — PnSmnPc Vm, n (z ICr 



(22) 



where p„ are non- negative c-numbers. Eq. (|22|) (or equiv- 
alently Eq.(pl|)) must be satisfied for some subset of ef- 
fects An (An) which form an operator-sum representation 
for the operation £. We emphasize that the An opera- 



tors in Eq.(|2l|) and the An operators in Eq.(22) are not 
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necessarily the same, and similar distinction holds for 
the reversible subsets. For simplicity of notation, we will 
drop the tilde even when we refer to the operators and 
reversible subset in Eq.(|2^). 

Eq.(p2[) is equivalent to the usual (exact) orthogonal- 
ity and non- deformation conditions for a non-degenerate 
code with logical states \ci), and these conditions can be 
stated concisely as: Vi, j. 



{Ci\Al^An\Cj) = S^jSmnPn n & K-r 



(23) 



where p„ are the non-negative c-numbers in Eq. (|T]) , and 
they represent the error detection probability of the effects 

An- 

When Eq.(p^ is satisfied, the effects An have polar 
decompositions 



AnPc = x/p^UnPc Vn e JCr 



(24) 



where the J7„'s are unitary. Note that PcUlUmPc = 
SnmPc is required for syndrome detection to be unam- 
biguous. This is simple to see by writing the recovery 
operation TZ as 



nip) = J2 ^kpRi + PepPe , 



(25) 



where Rk = PcUl is the appropriate reversal process for 
each effect Ak [k e /Cre), and Pe = I — J2keiCr 
When we apply TZ on £{\Tpin) {ipin\), where \Tpin) is a pure 
input state, the first term in Eq.(p5[) is proportional to 
|V'm)('0m|i aiid represents the recovered state. TZ{p) is 
trace preserving as a result of orthogonality of different 
UkPc- 

It is convenient to define F'''^* = J2neK P" 
total detection probability for the reversible subset. It is 
immediate that to achieve a desired fidelity J-', we have 
to include in Are a sufficient number of highly proba- 
ble effects so that P'^'^* > T . Thus, Are. is also a high 
probability subset. 

Now we generalize Eq.(|2^)-(|2^) based on the follow- 
ing assumption: the error is parameterized by a certain 
number of small quantities with physical origins such as 
the strength and duration of the coupling between the 
system and the environment. For simplicity, we consider 
only one-parameter processes, and let e be the small pa- 
rameter. For example, e can be the single qubit error 
probability. Suppose the aim is to find a code for a known 
E with fidelity: 



T>\- 0{e 



t+i\ 



(26) 



In the new criteria, it is still necessary that P'*'^* > !F, 
that is, Are has to include all effects An with maximum 
detection probability max|^.^)gc tr(|V'm)(V'™|^l^n) ~ 
C'(e^), s < t. However, it is not necessary to recover 
the exact input state; only a good overlap between the 



input and output states is needed. In terms of the condi- 
tion on the codeword space, it suffices for the effects to be 
approximately unitary and mutually orthogonal. These 
observations can be expressed as relaxed sufficient con- 
ditions for error correction. Suppose 



AnP 



Un 



PcAnAnPc ■ 



(27) 



is a polar decomposition for An- We define c-numbers 
Pn and A„ so that Pn and PnXn are the largest and the 
smallest eigenvalues of PcAl^AnPc, considered as an op- 
erator on C. The relaxed conditions for error correction 
are that: 



P„(l-A„) <0(e*+i) VnelCre 

PcUlUnPc = 5„,nPc ■ 



(28) 
(29) 



Note that when A„ = 1, Eq.(g7|) - (gOD reduce to 
the exact criteria. In the approximate case, P'''^* = 
Z^neK^e *^(l''^™)*^^™l"^n^") is not a constant, but de- 
pends on the input state IV'm)- Since Are includes enough 
effects so that P'^"* > 1 - 0(e*+i), when Eq.(|2|) is 
satisfied, we also have J^neic Pn ^ ^ — 0(e*^^) and 
E„e^..P«A„>l-0(e*+i). " 

We now prove that X^nex; Pn^n is a lower bound on 
the fidelity. Defining the residue operator 



= V PcAiAnPc - VPnKPc ■ 
we find, for the operator norm of 7r„, 
< |7r„| < VP^- \pPnKi-, 

and 

AnPc = Un{\/pnXnI + TTn)Pc ■ 



(30) 

(31) 
(32) 



The sufficiency of our conditions to obtain the desired 
fidelity may be proved as follows. Though Eq.(p^) is not 
satisfied, as long as Eq.(29) is true, we can still define the 



approximate recovery operation 

n{p) = RkpRl + PepPe 



(33) 



with Rk = PcUl as the approximate recovery operation 
for Ak and Pe defined as in the case of exact error cor- 
rection. For a pure input state \ipi7i){''Pin\, applying Tl{p) 
on f (|V'in)(V'in|)j and ignoring the last term which is pos- 
itive definite produces an output with fidelity 



= min tr(|-0„)(?/;m|7^(£(|V'm)(V'm|))) 

\{^^n\UlAn\^p^n)\'' ■ 



(34) 



Omitting all terms for which k ^ n and applying Eq. ( [32 
gives 
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T> min KV'mlVPnAn + 7r„|V'ii 



> 



(35) 



where in the last step, we have used Eq.(|3l|). Hence, the 
fidehty is at least X^neK Pn^n > 1 — 0(6*+^) and the 
desired fidelity order is achieved as claimed. 

An explicit procedure for performing this recovery is 
as follows. First, a measurement of the projectors Pk = 
UkPcUl is performed. Conditional on the result k of 
the measurement, the unitary operator Uk is applied to 
complete the recovery. 

Note that when the exact criteria hold, it is not neces- 
sary for the set A to form an operator sum representation 
for the error process; they may instead give a generalized 
description of the quantum process such as 



E 



(36) 



where Xmn is a matrix of c- numbers |15[ |. In this case, 
approximate criteria can be obtained straightforwardly 
along the lines we have described for the operator sum 
representation (where Xmn is diagonal). The main is- 
sue is to ensure that interference terms coming from off- 
diagonal terms in Xmn are sufficiently small. We will not 
describe this calculation here. Instead, we return to our 
four-bit code, and analyze it using the new criteria. 



IV. 4-BIT CODE REVISITED 

In terms of the approximate quantum error correction 
criteria Eqs.(|2^p9|), we may understand why our four-bit 
amplitude damping quantum code works as follows. Let 
the orthonormal basis for the Hilbert space be ordered 
as: 



10000), 10011), |1100), |1111), 10111), 10100), . 



(37) 



We include only the first few non-zero rows and columns 
in our matrices. The projection operator |0l)(0l| + 
|1l)(1l| onto the two dimensional codeword space C is 



Pr 



10 1 
110 
110 
10 1 



(38) 



The first effect operator (no loss of a quantum to the 
environment) is 



IQOOO 



10 

1-7 

1-7 

(1-7)^ 



(39) 



The eigenvalues of -Pc ^oooo^oooo-Pc are (1 — 7)^ and 
^(1 -I- (1 — 7)**). Interested readers can check for them- 
selves that 



^0000 

J/oooo \{1 - 1)1 + (7' + O(7*))^ooool Pc 



(40) 



(the order of 7 in ttoooo is factored out of ttoooo) with the 
choice: 



ttoooo — 



'''0000 — 





1 
1 



cos(f — 



sin(6l - f ) 

10 





1 



(41) 



(42) 



where tan 6* = (1 — 7)^. The exact criteria are not satis- 
fied, as ^0000 ^0000 has different eigenvalues. How- 
ever, the difference is of order 0(7^) and thus the relaxed 
condition Eq.(28) is satisfied. 

For the second effect (loss of one quantum from the ni 
mode), we have 



A 



1000 



(1-7)2/7 




















































































1-7 














1 












(43) 



The eigenvalues of Pc^iooo^iooo^b are 7(1 — 7) and 
7(1 — 7)^. The difference is (27^ — 7^)(1 — 7). We have 
the decomposition 



(1 - 7)7 



U 



1000 



(1 - 7)/ + 77riooo 



Pc 



(44) 



f^iooo — 



1 

1 

1 

1 

1 

1 



1 








1 











1 


1 

1^ 



























1 

V2 








1 

V2 




















1 




















1 



'''1000 


















1 

\/2 

















1 



(45) 



Other one loss cases are similar. The UkPc have non- 
zero entries in different rows, and are orthogonal to each 
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other. Hence, all the approximate code criteria are sat- 
isfied. Using these explicit matrices, we obtain 

7^(f(|V'^„)(^/'^n|))« ^ ^fe | V-™) (^/'^n 1 4 ^^fe^C 

- (1-372)|V',„)(V^,„| + ... . (46) 

The fidelity is thus at least 1 — 87^, and is of the desired 
order. 

The recovery procedure suggested in Eq. (^) contrasts 
with the decoding and recovery circuits in section 11 A , It 



is an interesting exercise to check that the composition of 
the operations in Fig. ^ and|^B, followed by re-encoding 
the recovered qubit has the same effect on C as applying 
C/qqqq for recovery. For the case in which an emission 
occurs in the first qubit, the composition of operations 
in Fig. 1^ and ||C, followed by re-encoding the recovered 
qubit has the same effect on C as preferentially damp 
out the jria) — |0) component followed by applying C/j'^qoo 
for recovery. Note that it costs 27^ in the fidelity for 
removing the distortion. 



JV, 



(50) 



That is, instead of requiring the individual number of ex- 
citations in each QCS to be balanced, it is sufficient to 
balance the average number of excitations over the QCS 
in each codeword. This provides an alternative expla- 
nation of why the known five-bit perfect quantum codes 

IS 

|0l) = 100000) + 111000) - llOOll) - 101111) 
+ 111010) + 100110) + 101101) + iioioi) 

|li) = 111111) - 100011) + 101100) - iioooo) 

- 100101) -I- 111010) + llOOlO) - 101010) (51) 

work for amplitude damping errors (as described in 
Eq.(^): although the codes do not satisfy the exact non- 
deformation criterion, Eq.(^) for the error representa- 
tion of Eqs.([^^, they do satisfy the approximate ones 
leading to Eq.(5Q). 



V. APPLICATIONS TO OTHER CODES 

Our approximate criteria may also be used to simplify 
code construction using a non-Pauli error description ba- 
sis. For example, consider the bosonic quantum codes 
for amplitude damping pO[ | (these are codes which utilize 
bosonic states |0), • • • \n) instead of qubits). For logical 
states |ci), |c2), ... of the form: 

\ci) = \/Mlkill»^12 • • ■ niyn) 
+ ^^\n2in22 ■ ■ ■ n2m) 

H 

+ ^/lWl\nNll'^Nl2 ■ ■ -nNim) , (47) 

the original non-deformation conditions for correcting up 
to one photon loss require the following to be constant 
for all logical states: 

(q|4Ao|q)=^(1-7)^^'M. (48) 

i=l 
Ni 

(q|4...i...o^o...i...o|q) = 5](1 - i)'''''-'ifi.n,, . (49) 

In the above, RSi — J2JLi ''^ij the total number of 
excitation in the i*'' quasi-classical state (QCS) in |q). 
It is difficult to find a solution for Eq.(|49|) when RSi is 
not constant for all i. That is the reason why previous 
work suggests using only QCS with the same number 
of excitations. 

With the new criteria, it suffices for the following to 
be constant for all logical states: 



VI. CONCLUSION 

We have shown by an example that choosing an ap- 
propriate error basis can potentially reduce the num- 
ber of qubits and other requirements in coding schemes. 
We also suggested a method to enable code construc- 
tion without the Pauli basis to be done more easily. We 
believe that much more can be done along these two di- 
rections. The new sufficient criteria for error correction 
are far from being necessary. Our results show that it is 
worthwhile to look for better approximate conditions or 
conditions which are both necessary and sufficient. Ap- 
proximate error correction is particularly interesting be- 
cause it is a property with no classical analogue, as it 
makes use of slight non-orthogonalities possible only be- 
tween quantum states. It would be especially useful to 
develop a general framework for constructing codes based 
on approximate conditions, similar to the group-theoretic 
framework now used to construct codes which satisfy the 
exact conditions. 
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